Ferromagnetism and Lattice Distortions in the Perovskite YTiOs 
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The thermodynamic properties of the ferromagnetic perovskite YTiOs are investigated by thermal 
expansion, magnetostriction, specific heat, and magnetization measurements. The low-temperature 
spin- wave contribution to the specific heat, as well as an Arrott plot of the magnetization in the 
vicinity of the Curie temperature Tc — 27 K, are consistent with a three-dimensional Heisenberg 
model of ferromagnetism. However, a magnetic contribution to the thermal expansion persists well 
above Tc, which contrasts with typical three-dimensional Heisenberg ferromagnets, as shown by 
a comparison with the corresponding model system EuS. The pressure dependences of Tc and of 
the spontaneous moment M s are extracted using thermodynamic relationships. They indicate that 
ferromagnetism is strengthened by uniaxial pressures p || a and is weakened by uniaxial pressures 
p || b, c and hydrostatic pressure. Our results show that the distortion along the a- and b- axes is 
further increased by the magnetic transition, confirming that ferromagnetism is favored by a large 
GdFe03-type distortion. The c-axis results however do not fit into this simple picture, which may 
be explained by an additional magnetoelastic effect, possibly related to a Jahn- Teller distortion. 



PACS numbers: 75.30.-m,75.50.Dd,75.50.Ee,75.80.+q 



1. INTRODUCTION 

ABO3 perovskites exhibit a large variety of electronic 
and magnetic properties 1 . The titanate family ATiOs re- 
cently attracted particular interest, since YTi03 orders 
ferromagnetically below the Curie temperature Tc — 27 
K, whereas LaTiOs orders antiferromagnetically below 
the Neel temperature T N ~ 150 K±&2£&&. In these 
systems, the 5=1/2 spins localized on the Tz 3+ ions 
are responsible for the magnetic properties. A change 
from ferromagnetism to antiferromagnetism can be con- 
tinuously tuned by varying the lanthanum concentration 
x in the alloys Yi-^La^TiOs, or by changing the lan- 
thanide A (A = Yb^La) in the undoped ATiOgiSa&l 
A GdFe03-type distortion is driven by ion-size mismatch 
and comprises rotations of the Ti06 octahedra. It is 
responsible for the distorted structure of the ATiOs cry- 
tals, with the space group Pbnm. This distortion is more 
pronounced in YTi03 than in LaTiOs, being favored by 
smaller A 3+ ions (A = Y,La) 4,5 . In YTi03, an additional 
elongation, by about 3 %, of the Ti06 octahedra is ob- 
served. This distortion has been ascribed to staggered 
ordering of the t2 g orbitals (Ti 3+ ions)^^. The switch 
from antiferromagnetism to ferromagnetism in the ATi03 
perovskites is probably controlled by the extreme sensi- 
tivity of the magnetic superexchange interactions to the 
distortions of the lattice 4 ^^. However, the mechanism 
driving this transition is still a matter of considerable 



debate^^^ii^iiii 1 ^^. For a proper description of the 
magnetic properties, it is thus crucial to carefully con- 
sider their dependence on the lattice distortion. 

In this article, we present a study of the thermody- 
namic properties of YTi03. Experimental details will 
be given in Section [2j In Section [3l the specific heat, 
thermal expansion, magnetization, and magnetostriction 
data measured with magnetic fields applied along the 
easy c-axis will be presented. In Section [H these results 
will be discussed and compared to the behavior expected 
within a three-dimensional (3D) Heisenberg ferromag- 
netic model 14 . As a specific example, we will show data 
on the typical 3D Heisenberg system EuS 15 i 16 . In Section 
[5j the relation between the distortion and the magnetic 
properties will be discussed in the light of our results. 
The dependence of the distortion on the A 3+ ionic sizes, 
on uniaxial pressures, and on the temperature will be 
considered. 



2. EXPERIMENTAL DETAILS 

Single crystals of YTi03 were prepared by the float- 
ing zone method using a four-mirror-type infrared image 
furnace from Crystal System Corporation. More details 
about the crystal growth are given in Ref. [I?]- Two 
samples have been investigated and the measurements 
presented here were obtained on the sample with the 
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sharpest transition at Tq. This sample was cut so that 
its faces are perpendicular to the a-, 6-, and c-axes, its 
dimensions at room temperature being equal to L\ ~ 2 
mm, L® ~ 4 mm, and L[? ~ 3 mm along a, 6, and c, 
respectively, with a mass of 116 mg. Thermal expansion 
and magnetostriction were measured using a home-made 
high-resolution capacitive dilatomete r 18 i 19 , with temper- 
ature and field sweep rates of 20 mK/s and 0.5 T/min, re- 
spectively. Three sets of measurements were performed, 
where the length Li was measured along the a-, 6-, and 
c-axes (i — a, 6, and c, respectively). Specific heat and 
magnetization were measured using a Physical Properties 
Measurement System and a Magnetic Properties Mea- 
surement System, respectively (Quantum Design). For 
all measurements, the magnetic field H was applied par- 
allel to the easy axis c. The thermal expansion of EuS 
was measured using a 8 mm long single crystal grown 
from the melt by K. Fischer at the Forschungszentrum 
Jiilich, as described elsewhere 20 . 



3. RESULTS 

3.1. Specific heat and thermal expansion 

In Fig. [TJ the specific heat C p of YTiC>3 is shown in 
a C p versus T plot. Ferromagnetic ordering is charac- 
terized by an anomaly at Tq = 26.8 ± 0.3 K, defined 
at the minimum of slope of C P (T). In Fig. [2] (a), the 
variation with T of the relative lengths ALi / Li is shown 
for i = a, 6, and c, ALi/ Li being fixed to zero at room 
temperature. The linear thermal expansion coefficients 
ai = (l/Li)dLi/dT, with i = a, 6, and c, are extracted 
from these data and are plotted in Fig. [2] (b). The vol- 
ume change AV/V = ^ i=a b c ALi/Li and the related 
volume thermal expansion coefficient oty = (l/V)dV/dT 
are also shown in Fig. [2] (a) and (b), respectively. As 
seen in Fig. O changes of slope in L a , L5, L c , and V are 
induced at Tc, leading to a positive anomaly in a a and 
to negative anomalies in 0^5, a c , and ay- The volume de- 
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FIG. 2: Variation with T (a) of the relative lengths ALi j Li 
for i = a, b, and c, and of the relative volume AV/V, (b) of 
the thermal expansion coefficients on for i = a, 6, c and V. 

crease below Tc, which is similar to the invar effect, will 
be further related to the negative hydrostatic pressure 
dependence of Tc. From the thermal expansion data, we 
extract a Curie temperature Tc = 26.8 ± 0.05 K at the 
extremum of slope of c^(T). This is in good agreement 
with prior observations 17 . At temperatures sufficiently 
higher than Tc, L a , L c , and V increase with T, while T5 
decreases with T (Fig. [2] (a)). This leads to the posi- 
tive values of a a , a c , and ay and to the negative values 
of ajj observed at high temperatures in Fig. [2] (b). The 
anisotropy of at high temperatures is a consequence of 
the lattice distortions, which will be discussed in Section 
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FIG. 1: Variation with T of the specific heat C p of YTi0 3 . 



3.2. Magnetization and magnetostriction 

In Fig. [3l the magnetization versus field M(H) is 
shown for H || c in a log- log plot, at T = 1.8 K and 
T = 26.5 K ~ T c . At T = 1.8 K, a linear increase of 
M(H) is obtained for fi H < fi H* ~ 0.06 T and is 
related to the alignment of ferromagnetic domains. For 
H > H*, the spins are aligned parallel to H and the 
magnetization M reaches M s ~ 0.8 Hb- In this regime, 
a slight increase of M(H) is observed. Indeed, M s is not 
yet fully saturated and is somewhat smaller than the full 
moment of 1 hb expected for the S = 1/2 Ti 3+ ion a 14 i 21 . 
At T ^ Tc, M increases first almost linearly with for 
fi H < fioH*' ~ 0.02 T, and then varies as M ex H 1 / 6 , 
with S = 4.47 ±0.2, for fj, H*' < fi H < 1 T. This power 
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FIG. 3: Field-dependence of the magnetization M, on a log- 
log scale, at T = 1.8 K and T — 26.5 K ~ T c . 
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FIG. 4: Field-dependence of the magnetostriction A^, on a 
log-log scale, for i = a, b, and c, at T = 2.5 K and Tc = 26.7 
K (X a is plotted with a minus sign). 



law, observed above 1 T, is characteristic of the critical 
ferromagnetic regime. Deviations are observed when M 
becomes close to M s . 

In Ref. [I?], low temperature magnetization measure- 
ments were reported on the same sample as in the present 
work and moments of 0.84 and 0.82 \±b were found for a 
magnetic field of 7 T applied along b and c, respectively. 
These data agree well with the saturated magnetic mo- 
ment of 0.84 ii b reported in Ref. [2l| andl22l and with the 
moments of 0.83 ± 0.05 and 0.84 ± 0.05 fi B , for the a- 
and c axis, respectively, determined by recent Magnetic 
Compton Profile experiments 23 . 

In Fig. [H the magnetostriction coefficients A^ = 
(l/Li)dLi/d(/jLoH) are plotted as a function of H on a 
log-log scale, at T = 2.5 K and Tc = 26.7 K, with i = a, 
6, and c, and H || c. For all temperatures and magnetic 
fields, A5 and A c are positive while A a is negative. For 
the three configurations at T = 2.5 K, \Xi(H)\ increases 
for fioH < fioH* ~ 0.06 T and is almost constant for 
H > iJ*, when the domains are aligned. At Tc, \Xi(H)\ 
increases before reaching a maximum at fioH*' 0.02 T. 
For fioH > fioH*', a critical regime is observed, where 
\Xi(H)\ ex if-V^ with SI = 2.64, 3.06, and 2.59 ± 0.1 
for z = a, 6, and c, respectively. While the power law is 
followed up to almost 10 T in A5, deviations are found 
for fioH > 1 T in A a and A c . 



4. FERROMAGNETIC PROPERTIES 
4.1. Low temperature spin waves 

In Fig. [5] (a), the specific heat of YTi03 is plotted 
in a log-log plot of C p /T versus T. Well below Tc, 
the phonon contribution can be neglected and the sig- 
nal, which varies as C P (T) oc T 1A up to 10 K, is be- 
lieved to be only magnetic. This power law is compatible 
with isotropic 3D Heisenberg ferromagnetic spin waves, 
for which a T 15 law would be expected, and is thus in 
good agreement with the spin-wave dispersion observed 
by neutron scattering 14 . The slight deviation from the 
T 1 - 5 law might result from a small spatial anisotropy of 
the exchange, the spin-wave contribution to the specific 
heat of a ferromagnet varying as C P (T) ex T d / 2 , where 
d is the dimensionality of the exchange. It may also 
be related to possible additional antiferromagnetic spin 
fluctuations and/or to spin anisotropics originating from 
spin-orbit coupling (see below). 

In Fig. 0(b), the thermal expansion is plotted in a log- 
log plot of \cti\ versus T, for i = a, 6, and c. Power laws 
oti oc T 1,9 are found up to almost 20 K for i = a and c, 
while no clear power law is observed for i = b. As simple 
3D Heisenberg ferromagnetic spin waves should lead to 
oti (x C p (x T 1,5 , the different T-dependences of C p and 
oti reported here may result from anisotropic exchange 
interactions, which lead to weak additional magnetic 
Bragg reflections due to canting of the ferromagnetic 
moments^. Weak low-energy spin fluctuations around 
these wave vectors may contribute to the deviation of the 
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FIG. 5: (a) Specific heat and (b) thermal expansion of YTiOs, 
in Cp/T and \m\ versus T on log-log scales, respectively. 



temperature dependence of the low-temperature specific 
heat and thermal expansion from the predictions of a sim- 
ple ferromagnetic Heisenberg model 24 . Further studies of 
the spin wave spectra (e.g. by neutron scattering) and 
detailed calculations are needed for a quantitative expla- 
nation of the results obtained here. Systematic studies by 
specific heat and thermal expansion of the ATi03 family, 
such as the work initiated in Ref. 0, may be of impor- 
tance to understand the evolution of the low temperature 
magnetic properties. 



mal expansion (Section [3+]) . From Equation (pQ): 

M(H : T C ) octf 1 /^ (2) 
with 5=(/? + 7)//3. (3) 

Using the exponents a and (3 obtained with Arrott's 
method, we calculate the critical exponent S = 4.76 + 0.2, 
which agrees favorably with 5 — 4.47 ± 0.2 directly ob- 
tained from the fit by a power law of M(H) at 26.5 K 
(Section [321). 

In Table H the exponents (3, 7, and S expected for dif- 
ferent classes of universality 27 are listed for comparison. 
The exponents extracted from the Arrott plot of the mag- 
netization of YTiC>3 are rather close to those of the 3D 
Heisenberg universality class. A similar plot was made 
in Ref. [28| using 3D Heisenberg exponents, but without a 
preliminary fit of the M(H, T) data as done here. As al- 
ready inferred from the behavior of the low temperature 
spin waves, the critical behavior of the magnetization, 
too, is thus consistent with a 3D Heisenberg picture of 
ferromagnetism (cf. Section RTT1 and Ref. Il4l ). 

In the following, the power laws reported for the mag- 
netostriction at Tc (Fig. 2j) are related to the critical 
power law of the magnetization (Fig. [3]) and to the crit- 
ical exponents a and (3. Using the Maxwell relation: 



A, 



1 dLi _ dM 
Li d(fjb H) dpi ' 



(4) 



the magnetostriction coefficients can be expressed as 
functions of the uniaxial pressure dependences of the 
magnetization. Assuming that dci/dpi = 0, the deriva- 
tive of Equation (pQ) leads to, at T = Tc: 



Xi (H,T c ) = -A^H-V*' 
dpi 



(5) 



with A 



\ ll8 'c 2l jb and <y / =(/? + 7 )/(l-/3). (6) 



From the exponents a and (3 obtained by the Arrott fit 
of the magnetization, a critical exponent S f = 3.07 ± 0.5 



4.2. Arrott plot and critical fluctuations 

To analyze the critical ferromagnetic regime, the mag- 
netization was measured as M versus H at several tem- 
peratures close to Tc. In Fig. [6l an Arrott plot of these 
data is shown as M 1 ^ versus (H/M) 1 ^, for 25.9 < T < 
27.1 K. The critical exponents (3 = 0.392 ± 0.05 and 
7 = 1.475 ±0.1 used in this plot were determined from a 
fit of M(H,T) using Arrott's equation of state 2 ^^: 

M 1 ^ = Cl - pa (T-r c ), (1) 




(jU H/M) ((Oe.mol/emu) 



for 0.02 < fi H < 1.3 T and 25.9 < T < 27.1 K. The 
value of Tc = 26.44 + 0.05 K obtained from this fit agrees 
well with that obtained from the specific heat and ther- 



FIG. 6: Best Arrott plot of the data, with (3 = 0.392, 7 = 
1.475, and Tc — 26.44 K. The dotted line indicates the critical 
regime at Tc, associated with the exponent 5 = 4.76. 
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TABLE I: Critical exponents /3, 7, 8, and 8' obtained here for 
YTi03 and expected for different universality classes^. 

Critical exponents f3 7 5 8' 

YTi0 3 (best fit) 0.392 (50) 1.475 (100) 4.76 (20) 3.07 (50) 
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is estimated by the blue area. The estimates of the mag- 
netic contribution to the thermal expansion of YTiC>3 
and EuS are plotted in Fig. [71 (c). In this plot, the mag- 
netic thermal expansion coefficient a mag is normalized by 
its maximal value ano ^ the temperature T is nor- 

malized by Tc- Fig. 0(c) indicates that, in YTiC>3, the 
magnetic signal has a significant weight above Tc and 
extends up to about 5 x Tc while, in EuS, it has most 
of its intensity below Tc and vanishes completely above 
about 2 x Tc. Thus, the magnetic fluctuations of YTiC>3 
cannot be described as those of a simple 3D Heisenberg 



is expected to characterize the magnetostriction at Tc- 
This value is in good agreement with the values 2.59 < 
5^ < 3.06 determined from the fits of Xi(H) at Tc (see 
Section [3^2]) . For each set of (/3/y), the corresponding 8' 
values are also given in Table H These values are compat- 
ible with a 3D Heisenberg scenario of ferromagnetism for 
YTiOs associated with 8' = 2.77. The slight variations 
of 8[ with i are not understood and may result from var- 
ious secondary effects (anisotropic energy scales, defects 
etc.). 

To our knowledge, YTi03 is the second ferromagnetic 
system known, after the itinerant ferromagnet Ulr— , 
where a critical power law is reported in the magne- 
tostriction at Tc. We believe that such an effect is quite 
general and should be present in most ferromagnet s, once 
the field and temperature ranges are properly chosen. We 
note that, more than 60 years ago, Belov has theoretically 
predicted a similar law 30 , but only within a mean- field 
approach, which corresponds to (3 = 0.5, 7 = 1, and 
8 = 8' = 3 (cf. Tablelj). Our approach is more gen- 
eral and permits to obtain the critical exponent 8' for 
each combination of (/3/y) and thus, for each universality 
class. 



4.3. High-temperature magnetic signal: deviation 
from a pure 3D Heisenberg ferromagnet 

The distortion of the lattice induces the strong 
anisotropy observed in the high-temperature thermal ex- 
pansion of YTi03 (see Fig. [2]). This anisotropic lattice 
signal must be properly taken into account to extract the 
magnetic contribution to the thermal expansion. Since 
a a and a c are rather close above 100 K, we assume that 
the lattice contributions to a a and a c are similar, so that 
Aa = a a — a c , which is plotted in Fig. [71 (a), can be con- 
sidered as a signal representative of the magnetic thermal 
expansion. In this plot, the thin red line is a guide to the 
eye indicating the "remaining" non-magnetic background 
and the yellow area corresponds to the estimated mag- 
netic contribution. In Fig. 0(b), the thermal expansion 
a of cubic EuS 31 , which is known as a prototype of 3D 
Heisenberg ferromagnetis m 32 i 33 i 34 , is shown for compari- 
son. The thin blue line is a guide to the eye indicating the 
non-magnetic background and the magnetic contribution 



YTiOa 




T/T c 

FIG. 7: Temperature dependence of (a) (Aa = a a — a c ) for 
YTiOs, and of (b) a for EuS. In both plots, the thin line is a 
guide to the eye indicating the non-magnetic background, the 
colored area is the magnetic contribution deduced from the 
raw data after substraction of this background, and the ver- 
tical dotted line indicates the ferromagnetic transition tem- 
perature, (c) Magnetic contribution to the thermal expansion 
estimated for YTiOs (red line) and for EuS (blue line), in a 
normalized Qt ma9 /ot™ 9 x vs T/T c plot. 
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FIG. 8: Temperature dependence of the magnetic contribu- 
tion C™ a9 to the specific heat, in a C™ a9 /T versus T plot. 
The inset shows the scaling of C P (T) (red line) and av(T) 
(green line), which permitted to estimate the non-magnetic 
contribution (blue line). 



ferromagnet. 

As shown in the inset of Fig. [51 the specific heat 
C P (T) and the volume thermal expansion ay{T) can be 
scaled at high temperatures using an empirical parameter 
y = 3.96 x 10 -6 mJ/mol defined by C p = y x ay. Assum- 
ing that the magnetic contribution to C p and ay is neg- 
ligible above 200 K, and that there is a single Gruneisen 
parameter associated with the phonons for the [0,300 K] 
range, we can estimate the phonon contribution to the 
specific heat by C% h = C p — z x (C p — y x ay), z = 0.2 
being adjusted so that no anomaly remains at Tc- The 
main frame of Fig. [8] shows the resulting estimate of 
the magnetic contribution C™ a ^(T) to the specific heat. 
This plot confirms the conclusions from Fig. [3 i.e. that 
a magnetic contribution is present up to more than 100 
K. Integration of the estimated magnetic heat capacity 
leads to the magnetic entropy AS mag ~ 4.5 J/mol.K. 
This entropy is roughly equal to the full spin entropy 
A5™f = R\n2 ~ 5.8 J/mol.K expected for S = 1/2 spin 
system; the fact that AS ma9 is about 20 % smaller than 
A5^f may be imputed to the experimental error. 

In Fig. [51 the estimated magnetic contributions to 
the thermal expansion of mafir (T) and to the specific heat 
C™ a9 (T) are plotted together; above 50 K, the different 
shapes of the signals indicate the limit of the methods 
used here. Both plots indicate clearly the presence of 
a magnetic signal at temperatures well above ferromag- 
netic ordering. The origin of this behavior, which is not 
expected for usual 3D-Heisenberg ferromagnets (see Fig. 
[7]), is not yet understood. In principle, a modified con- 
ventional spin-only fluctuation model, for example with 
competing (and possibly low-dimensional) antiferromag- 
netic and ferromagnetic interactions, could describe this 
anomalously high-temperature magnetic signal. How- 
ever, this is in apparent contradiction with the magnon 
spectra reported by neutron scattering, which do not ex- 
hibit pronounced deviations from the predictions of a 3D 



Heisenberg model with nearest-neighbor interactions^. 
An alternative explanation of the extended magnetic fluc- 
tuation regime could be offered by spin-orbital fluctua- 
tions models, where an energy scale significantly exceed- 
ing the magnon bandwidth (in Ref. [4HI42L orbital fluc- 
tuations were associated with an excitation at about 250 
meV) could actuate ferromagnetic fluctuations at tem- 
peratures between Tc and room temperature^. Further 
work is required to ascertain whether a quantitatively 
consistent picture of the spirit and orbita l 41 ! 42 excitation 
spectra and thermodynamics of YTi03 can be obtained. 



5. COUPLING BETWEEN THE LATTICE AND 
THE MAGNETIC PROPERTIES 

5.1. Uniaxial pressure dependences - Comparison 
with LaTiOs 

Fig. [9] shows the anomalies at Tc in the specific heat 
and in the thermal expansion of YTiOs. These anoma- 
lies are typical of a second-order phase transition, whose 
jumps are estimated as AC P = 2.3 ± 0.1 J/molK in the 
specific heat and as Aa a = 6.6 ±0.3 x 10 -6 K _1 , Aa b = 
-3.4±0.2 x 10~ 6 K-\ and Aa c = -4.7±0.2 x 10~ 6 K" 1 
in the thermal expansion. Using the Ehrenfest relation: 

dTc = AajVmTc () 
d Pi AC P ' U 

where V m = 3.46 x 10 -5 m 3 / mol is the molar volume and 
Pi a uniaxial pressure applied along i (i = a, 6, and c), 
we extract the uniaxial pressure dependences dTc /dpi 
reported in Table [Til The sum of the three uniaxial pres- 
sure dependences of Tc gives the hydrostatic pressure 
dependence dT c /dp h = -6.0 ± 0.6 x 10~ 2 K/kbar. As- 
suming that, in YTiC>3, dTc/dph remains constant under 
hydrostatic pressure, ferromagnet ism may be destroyed 
above p c ~ 400 kbar. 

For comparison, the uniaxial pressure dependences of 
Tat, for the antiferromagnet LaTiOs, are also listed in 
Table [III To calculate them, Tc in Equation (0) was 
substituted by T N ~ 146 K, AC P = 10 J/molK, Aa a = 
-5 ± 0.5 x 10~ 5 K-\ Aa b = 5 ± 0.5 x 10~ 5 K _1 , and 



TABLE II: Uniaxial and hydrostatic pressure dependences of 
Tc and M s for YTi03 and of Tn for LaTiOs. The ratio pi of 
the pressures dependences of Tc and M s is given for YTiOs- 
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FIG. 9: Ferromagnetic ordering anomaly (a) in the specific 
heat and (b-d) in the thermal expansivity coefficients along 
a, 6, and c. In these curves, the jumps AC P , Aa a , Aa^, and 
Aa c at the ferromagnetic ordering are indicated by arrows. 

Aa c ~ being estimated from Ref. |4C|. W hile dTjsf/dp c 
is too small to be extracted from Ref. [4(J, dT^/dpa and 
dT^/Opb are such that dTjy /dp a = —dT^/dph < 0. For 
p || a, b, the uniaxial pressure dependences of Tc and 
Tat, for YTi03 and LaTiOs, respectively, have thus op- 
posite sign. The effects of pressure along c are such that 
dTc/dpc < and dT N /dp c ~ 0. In Sections O and 
15. 31 the uniaxial pressure-dependences of Tc and T/v in 
YTi03 and LaTi03, respectively, will be interpreted as 
resulting from pressure-induced modifications of the dis- 
tortion. 

Well below T c and for fi H > /i #* = 0.06 T, 
M ~ M s 0.8 \±b (Fig. [3]) and the magnetostriction 
coefficients of YTiOs are almost constant, having the val- 
ues A a ~ -3.5 ± 0.5 x 10- 6 T- 1 , X b ~ 3.5 ± 0.5 x lO" 6 
T- 1 , and A c ~ 1.0 ± 0.2 x lO" 6 T" 1 (Fig. g). Using the 
Maxwell relation given in Eq. (jlj), we extract: 

^ = -A,(T = 2.5K). (8) 
dpi 

The values of dM s /dpi, for i = a, and c, as well as their 
sum, the hydrostatic pressure dependence dM s /dph, are 
summarized in Table ITU 

For each i = a, 6, c, and h (h ^ hydrostatic), the pres- 
sure dependences dTc/dpi and dM s /dpi have always the 
same sign, being both positive for i = a and both nega- 
tive for i = b,c,h (Table [TTJ) . Ferromagnetic order is thus 
stabilized by uniaxial pressure p || a and is destabilized 
by uniaxial pressure p || b, c and by hydrostatic pressure. 
Consequently, the ratio p^, defined by: 



M s dTc/dp, 



(9) 



rz TcdMs/dp^ 

is always positive. As shown in Table [III we fi n d that 
pi is strongly anisotropic, being bigger when the z-axis is 
easier (c = easy, a = intermediate, and b = har d 17 i 21 ). 



Although YTiOs is a localized ferromagnet, its satu- 
rated moment, at about 5 T, is only 80 % of the fully 
saturated moment M^ ul1 — 1 \±b- A small canted anti- 
ferromagnetic moment Maf ~ 0.1 hb was reported by 
neutron scattering in YTiOa^ and explains partly why 
M s is reduced. In addition, the reduction of M s may indi- 
cate an enhanced phase space for quantum magnetic fluc- 
tuations. The question is whether this reduction comes 
from usual spin-only fluctuations, or if it results from 
more complicated fluctuations involving orbital degrees 
of freedomi&ii xhe high values of dM s /dpi in YTi0 3 
may result from the combination of two effects, which can 
be summarized as the uniaxial pressure-induced transfers 
of weight (i) between the ferromagnetic moment M s and 
the antiferromagnetic moment Maf and, (ii) between 
M s and some quantum magnetic fluctuations 5M. The 
second effect is similar to what happens in itinerant fer- 
romagnets, where M s is reduced by quantum fluctuations 
SM of the magnetic moment, and where the strong pres- 
sure dependences of M s are related to those of SM. 

A similar analysis as the one presented here was re- 
ported for the itinerant weak ferromagnet Ulr (Ref. I29I ), 
in the framework of the Moriya's spin fluctuation theory 
of itinerant magnetis m 38 ! 39 . By analogy, a spin fluctu- 
ation theory, adapted to the particular case of YTiOs, 
may be appropriate. 



5.2. Coupling between the GdFeC>3-type distortion 
and the magnetic ordering 

In the ATi03 perovskites, the GdFe03-type distor- 
tion comprises a combination of tilts and rotations of 
the Ti06 octaedra. This results in an orthorhombic 
structure, where b > aoV2 > c/y/2 > a, ao being the 
lattice parameter of an undistorted cubic structur e 43 ! 46 . 
In the alloys Yi^La^TiOa, La-substitution induces a 
decrease of the GdFe03-type distortion, which is be- 
lieved to control the change from ferromagnetism to 
antiferromagnetismi^i 3 -^ (Fig. [TO]) . This picture, in 
which ferromagnetism is favored by a large GdFe03-type 
distortion, is qualitatively confirmed by the increase of 
the distortion of the (a, b) plane induced below Tc (see 
Fig. [2]). In the following, we will further show that, for 
p || a, b, the uniaxial pressure dependences of Tc and 
T/v are mainly controlled by those of the GdFeOs-type 
distortion. 

In YTiOs and LaTiOs, the negative sign of dTc/dpb 
and the positive sign of dTjy/dpb (see Table HI]) , respec- 
tively, imply that a uniaxial pressure p || b can be seen 
as equivalent to La-doping (cf. the corresponding arrow 
in Fig. [TO]) . Conversely, the fact that dTc/dp a is positive 
while dT^/dpa is negative (Table HI]) implies that p || a is 
equivalent to Y-doping (see Fig. [10]). As p || b induces a 
compression along b and, because of elasticity, extensions 
along a and c, its effects are very similar to those of reduc- 
ing the GdFe03-type distortion. By analogy, p || a leads 
to a compression along a and to small extensions along b 
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GdFe0 3 -type distortion 



5.3. Distortion of the TiC>6 octahedra 



Y,.,La,Ti0 3 






/ AF Order 


T c 1 




FM Order\/ 





YTi0 3 


LaTiO 


p II a 






> 



p II b 

FIG. 10: Schematic magnetic phase diagram of Yi_ x La x Ti03 
(AF = antiferromagnetic and FM = ferromagnetic). The ar- 
rows indicate the effects of increasing the GdFe03-type dis- 
tortion and of applying uniaxial pressures p || a, b. 



and c, which is similar to increasing the GdFe03-type dis- 
tortion. Thus, we conclude that p || a and p || b induce 
an increase and a decrease of the GdFe03-type distor- 
tion, respectively, which are responsible for the various 
signs of dTc y N /dpi, for i = a, b. 

The uniaxial pressures dependences of Tc and XJv are a 
consequence of the high sensitivity of the superexchange 
interactions to the bond angles between the ions, whose 
positions are very sensitive to the pressure-induced modi- 
fications of the GdFe03-type distortion. The application 
of uniaxial pressures p || a, b, as well as the variation 
of the A 3+ ion, permits thus to tune the competition 
between the ferromagnetic and the antiferromagnetic ex- 
change interactions, via a change of the GdFe03-type 
distortion. 

However, a pressure-driven change of the GdFe03- 
type distortion cannot explain the results obtained for 
p || c, i.e. dTc/dp c < and dT N /dp c ~ 0. In- 
deed, dTc/dp c > and dT^/dpc < would be ex- 
pected if p || c merely modified the GdFe03-type dis- 
tortion (since p || c contracts c, it should increase the 
GdFe03-type distortion). Another mechanism, in addi- 
tion to the GdFe03-type distortion, is needed to under- 
stand the pressure dependences of Tc and T/v for p || c. 
In the next Section, we will show that a higher sensitiv- 
ity of the c-axis length to the intrinsic elongations of the 
octahedra may be the origin of this behavior. 



5.3.1. Microscopic description 

The lattice structure of YTiC>3 is represented schemat- 
ically in Fig. [Ill where the alternation of tilts and ro- 
tations of the Ti06 octahedra (in blue) is due to the 
GdFeC>3-type distortion. An additional distortion con- 
sists of an elongation of each octahedron along a partic- 
ular axis (orange arrows in Fig. [11]), and of contractions 
perpendicularly to this axis. In Ref. 45 ; 7 ; 8, the fact 
that the elongated axes vary from one site to another 
was ascribed to a staggered ordering of the ti g orbitals 
(Ti 3+ ions) via a collective Jahn- Teller effect. By ge- 
ometrical considerations, we can qualitatively estimate 
the macroscopic distortion induced by the elongations of 
the octahedra. As seen in Fig. [Ill the elongated axes are 
almost contained within the (a, b) plane, i.e. perpendic- 
ularly to the c-axis. This implies that the elongations of 
the octahedra induce a contraction of the c-axis. In the 
(a, b) plane, the elongated axes of two adjacent octahedra 
subtend an angle of about 60 °, so that the elongations 
and contractions of the different octahedra almost cancel 
each other. Since the projections of the elongated axes 
are larger along b than along a (the elongated axes sub- 
tend an angle of about 30 ° with b and of about 60 ° with 
a), we finally conclude that the elongations of the octa- 
hedra are responsible for a small elongation of b and for 




FIG. 11: Schematics of the lattice structure of YTiOs. The 
Ti 3+ ions are represented by grey spheres, the Y 3+ ions by 
green spheres, and the O 2- ions by red spheres. The Ti06 
octahedra are colored in blue, and orange arrows show their 
elongated direction, possibly due to the Jahn- Teller distor- 
tion. 
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a tiny compression of a, in addition to the main effect, a 
compression along c. 

5.3.2. Lattice parameters: comparison of the families 
ATi0 3 and AFe0 3 

Here we propose a method, based on a comparison of 
the lattice parameters a, b and c of the families ATiOs 
and AFeOs, to confirm the description made in Section 
15.3.11 of the effects of the elongation of the octahedra on 
a, b and c. Assuming that these elongations are related to 
a Jahn- Teller distortion 44 , the comparison of the lattice 
parameters of ATiOs and AFeC>3 can be justified by the 
fact that, contrary to Ti 3+ , Fe 3+ is not Jahn- Teller active 
so that AFeOs can be considered as a non- Jahn- Teller 
reference for ATi03. 

In Fig. [12] (a) and (b), the unit cell volume V and the 
lattice parameters a, 6, and c/y/2 are plotted versus the 



ionic radius of the A 3+ ions, for several compounds of the 
families ATi0 3 and AFe0 3 (A=Lu^La^^^^. The 
ionic radii of the A 3+ ions are taken from Ref. 0, assum- 
ing a number of 8 nearest neig hbors 6 . While LaTi0 3 and 
LaFe03 are almost undistorted (a ~ b ~ c/y/2), Fig. [12] 
(b) shows a strong distortion of the pseudo-cubic lattice 
in ATiOs and AFeOs, once A 3+ is smaller than La 3+ . 

In Fig. H21 (c), the unit cell volumes V ATi °3 Q f the 
ATi03 compounds are scaled empirically with the unit 
cell volumes y AFe °3 Q f the AFe03 compounds, using a 
scaling factor / = 1.01 defined by V ATi °3 = yAFe0 3 ^ 
f 3 . In Fig. [12] (d), the lattice parameters of ATi03 are 
scaled to those of AFe03 using the factor 1//. As the 
undistorted limit in the ABO3 perovskites corresponds 
to a cubic lattice parameter a^ B ° 3 = 2(ro + ^b), where 
ro and re are the ionic radii of the O 2- and B 3+ ions, 
respectively, we associate the empirical scaling factor / = 
1.01 to the ratio a^™ 3 /^ Fe ° 3 = 1.013, calculated with 
r = 1.35, r Ti = 0.67, and r Fe = 0.645 A 49 . 
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FIG. 12: Variations, for the ATiOs and AFeOs systems, of (a) the unit cell volume V and (b) the lattice parameters a, b, and 
c/y/2, as a function of the ionic radius of A 3+ . In (c) and (d), the volume and the lattice parameters of ATiOs and AFeOs are 
scaled together using the empiric factor / = 1.01. 



Since the ionic radii of Ti 3+ and Fe 3+ are very close, 
we assume that, for the two families, the GdFe03-type 
distortion induces similar variations of their lattice 
parameters in the scaled plot of Fig. [12] (d). Conse- 
quently, the elongations of the octahedra, which can 



be neglected in the non- Jahn- Teller compound AFeOs, 
might be responsible for the slight differences, in Fig. 
[T2l (d), between the scaled lattice parameters of the 
two families. This implies that the elongations of the 
Ti06 octahedra in YTiOa induce a decrease of c/y/2 by 
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about 0.5 A, accompanied by a smaller increase of 6, by 
about 0.3 A, and by no noticeable change of a. These 
conclusions, obtained using the scaled plot of Fig. [12] 
(d), confirm those deduced from geometrical arguments 
in Section 15.3. li 



5.3.3. Anomalous character of the c-axis? 

When T is reduced, the decrease of c is slowing down as 
the ferromagnetic transition at Tc is approached, which 
ends by an upturn below Tc, where c increases with de- 
creasing T (see Fig. [2] (a)). The behavior of the c-axis 
contrasts with those of the a- and b- axes, whose varia- 
tions are monotonic for 4 < T < 300 K and are ampli- 
fied below Tc (see Fig. [2] (a)). The distortion of the 
(a, b) plane results mainly from the GdFeC>3-type distor- 
tion, whose modifications also control the uniaxial pres- 
sure dependences of Tc and T/v for p || a, b (see Sec- 
tion 15. 2p . In Sections 15.3.11 and I5.3.2[ c was shown to 
be more sensitive than a and b to the elongations of 
the TiC>6 octahedra, possibly related to a Jahn- Teller ef- 
fect. The sensitivity of c to the distortions of the octahe- 
dra may be related to the anomalous uniaxial pressure- 
dependences of Tc and T/v for p || c (see Section [52]) , 
but also to the anomalous behavior of the c-axis in the 
spectral weight transfers of the optical conductivity^. 
Our findings are in apparent contradiction to theories 
according to which the Jahn- Teller distortion is an essen- 
tial prerequisite of ferromagnetism in YTiOs (Refs. 01!) • 
Rather, the properties of YTiC>3 seem to be reminiscent 
of those of Lay/gSrx/gMnOs, where a Jahn- Teller distor- 
tion is fully suppressed at the onset of ferromagnetism^. 

5.4. High-temperature extrapolation 



^ 










YTi0 3 







1000 2000 3000 4000 

T(K) 

FIG. 13: High temperature linear extrapolation of the lattice 
parameters a, 6, and c/y/2 of YTiOs. 

Equivalently to uniaxial pressures and variations of the 



A 3+ ion size, increasing the temperature leads to a reduc- 
tion of the distortion in the ABO3 perovskites. Indeed, 
the strong anisotropy of the thermal expansion coeffi- 
cients ai, shown in Fig. [2] (b), is governed by the effects 
of temperature on the crystal distortion. Assuming con- 
stant thermal expansion coefficients above room temper- 
ature, Fig. [13] shows high-temperature extrapolations, 
up to 4000 K, of the lattice parameters a, 6, and c/V^ 51 . 
This plot indicates that, in YTiOs, a cubic structure with 
a = b = c/y^ 43 might be recovered around 3000-4000 
K. However, this temperature scale, characteristic of the 
lattice distortion, is inaccessible since it is far above the 
melting temperature of YTiOs. 



6. CONCLUSION 

The thermodynamic study of the perovskite system 
YTiOs presented here allowed us to extract information 
about the ferromagnetic ordering and its coupling to the 
lattice distortions. While the low- temperature specific- 
heat data, as well as an Arrott plot of the magnetiza- 
tion close to Tc, are consistent with a 3D Heisenberg 
picture of ferromagnetism, deviations from this simple 
picture were observed in the thermal expansion data 
at low temperature, where an unexpected power law is 
found. Above Tc, a magnetic signal persists up to the 
remarkably high temperature of 5 x Tc- Further work 
is required to show whether models incorporating com- 
bined spin-orbital fluctuations, instead of spin-only fluc- 
tuations, could quantitatively describe this extended fluc- 
tuation regime. 

Ehrenfest and Maxwell relations enabled us to extract 
the uniaxial pressure dependences of the Curie temper- 
ature Tc and of the spontaneous moment M s , which 
indicates that ferromagnetism is stabilized by uniaxial 
pressures p || a and is destabilized by uniaxial pressures 
p || b, c and by hydrostatic pressure. We interpreted the 
uniaxial pressure dependences of Tc and M s obtained for 
p || a, b as resulting from uniaxial pressure-induced mod- 
ifications of the GdFeOs-type distortion. A high sensi- 
tivity of the c-axis to an additional distortion of the Ti06 
octahedra, possibly related to a Jahn- Teller effect, is be- 
lieved to be responsible for the anomalous uniaxial pres- 
sure dependences of Tc and M s observed for p || c. This 
confirms that both kinds of distortion play an impor- 
tant role for the formation of ferromagnetism in YTiOs, 
Jahn- Teller distortion being not a necessary condition 
for ferromagnetism in YTiOs. While the a- and b- axes 
are more sensitive to the GdFe03-type distortion, the c- 
axis is more sensitive to the elongations of the octahedra. 
Finally, a high-temperature extrapolation of the lattice 
parameters led to the onset of the distortion at a virtual 
temperature of about 3000-4000 K. These results might 
be considered to further develop models for the electronic 
properties of the titanates. 
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